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Abstract 



We observe critical phenomena in spherical collapse of radiation fluid. 
A sequence of spacetimes S[rj] is numerically computed, containing models 
(rj <C 1) that adiabatically disperse and models (i] 3> 1) that form a black hole. 
Near the critical point (rj c ), evolutions develop a self-similar region within 
which collapse is balanced by a strong, inward-moving rarefaction wave that 
holds m(r)/r constant as a function of a self-similar coordinate £. The self- 
similar solution is known and we show near-critical evolutions asymptotically 
approaching it. A critical exponent (3 ~ 0.36 is found for supercritical (ij > i] c ) 
models. 
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New solutions of Einstein's equations at the threshold of black hole formation have re- 
cently been discovered. These solutions exhibit nonlinear dynamical behavior closely anal- 
ogous to critical phenomena, and include power-law behavior, discrete scaling relations and 
a form of universality. All of these features were discovered in Choptuik's |I| studies of 
spherical collapse of massless scalar field. Abrahams and Evans found strikingly similar 
behavior in axisymmetric collapse of gravitational waves, indicating the behavior is a generic 
feature of gravity. Confirmation has recently been made |§ of some of Choptuik's results. 
In this Letter, we report observation of critical phenomena in spherical collapse of radiation 
fluid and show the asymptotic approach of near-critical spacetimes to a self-similar solution 
near the center of collapse. This new result is fundamental in two ways; the self-similar 
solution has been separately and precisely calculated [|J , and it exhibits local self-similarity. 

The process for finding the threshold of black hole formation and critical phenomena in 
gravitational collapse has been described elsewhere |l],|2],|5]|| . Adopting a physical model with 
metric g, matter fields 4>a, stress-energy T(c/)a), and symmetries, the future development is 
computed of elements of a family of Cauchy data, distinguished by a dimensionless strength- 
parameter 7] and a length (or mass) scale r . Subcritical models (77 1) have smooth future 
development, dispersing all mass-energy to infinity, while a black hole appears in supercritical 
models (77 ^> 1) with only some mass-energy escaping. As 77 is tuned in a search for the 
onset (at r/ c ) of black hole formation, black hole mass (for 77 > rj c ) diminishes as a power law: 
M B h = K{r] — rjc) 13 , with (3 > 0. In both scalar field collapse |l| and vacuum collapse [@J|, 
simulations give (3 ~ 0.37, results that appear to be universal (i.e., independent of the 
parameterization of the initial data). 

By implication, tuning the parameter to just the right value (77 = i] c ) creates arbitrarily 
small black holes and regions of arbitrarily high spacetime curvature from smooth Cauchy 
data. Hence, it is argued |J that the threshold of black hole formation is also the threshold 
of naked singularity formation. Spacetimes with rj = rj c are termed precisely critical. 

In any physical model, many families of Cauchy data may exist with a critical point and 
a precisely-critical spacetime. Each family of data will have a characteristic proper length 
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scale tq. Related to this scale is the limiting interval, lim \ T\ = \Tq\, as 77 — 7] c — > + , of proper 
time T of the observer O centered in the collapse. Let T = — T correspond to the initial 
data. Further, let r be a proper spatial distance from O. The central region of collapse 1Z 
(in spacetime) exists for r <C ro, \T\ <C \Tq\. Different precisely-critical spacetimes can be 
normalized, using scale invariance, to a common length scale r . These spacetimes will still 
differ on scales T ~ —To and r ~ ro but approach a (physical-model-specific but otherwise 
unique) self-similar solution within 1Z. Within 71, all knowledge of the initial data is lost up 
to a single, dimensional constant. 

Previous critical systems have displayed discrete self-similarities. Let q denote dimension 
and let length have q = 1. Allow the metric g to carry dimension q = 2 [[/J. For rj = i] c , within 
1Z a map will exist under which g — > g' = e~ 2A g, for a fixed value of A. These are echoes 
resulting purely from nonlinearities. The scale factors A are pure numbers that emerge from 
the dynamics, with e A ~ 30 in scalar wave collapse [Q and e A ~ 1.8 in gravitational wave 
collapse |@||. 

An analytic explanation of these phenomena would be important. This has been sought 
for scalar field collapse but the discrete self-similarity proved to be an impediment. The effort 
made obvious, though, that a local self-similarity would be more powerful, and suggested 
examining another system: perfect fluid collapse. 

Imagine a spherical configuration of radiation fluid confined within areal radius r . Let 
the total gravitational mass be M. Define a dimensionless control parameter 77 = 2M/ r [|J . 
When rj 1, pressure dominates and should cause adiabatic expansion and dispersal of the 
fluid. When 77 > 1, gravity should overwhelm pressure and spur the formation of a black 
hole. The onset of black hole formation should occur at some r\ c ~ 1. One might anticipate, 
with a fluid, that any self-similarity emerging near r] c will be local. 

A successful attempt was made to find a self-similar solution representing the threshold 
of black hole formation in perfect fluid collapse 0. Details of this solution will be given 
elsewhere; here we only describe the self-similar ansatz and a few numerical results. We 
assume a perfect fluid T^ v = pU^lf" +pg fiu , with total energy density p, isotropic pressure 
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p and four- velocity If 1 satisfying U^W 1 = — 1. A relativistic equation of state p = (7 — l)p 
is adopted. In this paper, we specialize to radiation fluid p — |p. Adopting (dynamical) 
Schwarzschild coordinates, the line element is 

ds 2 = -a 2 dt 2 + a 2 dr 2 + r 2 dVl 2 , (1) 

with a(r,t) the lapse function and a(r,t) the radial metric function. Constant proper 
sound speed c s = \Jpf~p = V7 ~ 1 suggests the self-similar ansatz: a(£), U r (£), p(r,t) = 
Q(^)/4nr 2 , a(r,t) = nrN(!;)/t, with self-similar coordinate £ = r/C(—t) n . Geometrically, 
w = (l/n)td/dt + rd/dr is the infinitesimal generator of the homothetic motion: <£ w g = 2g, 
£ W T = 0. The similarity exponent n is initially unknown, making the problem (in these co- 
ordinates) a type-2 self-similarity ||. With this ansatz, the fluid and gravitational equations 
reduce to ordinary differential equations. These are solved by demanding regularity at £ = 
(r = 0, t < 0), £ = 00 (r > 0, t = 0), and along the limiting, ingoing acoustic characteristic 
(sonic point) at £ = 1, across which the fluid might otherwise be discontinuous. The system 
is over-determined unless n is taken as an eigenvalue of the nonlinear system. An analogous 
type-2 similarity problem is Guderley's shock implosion problem [ |T0|j . 

The similarity exponent was found to be n ~ 1.1485 (7 = 4/3). The solution is depicted 
in Fig. 1. A physical singularity exists at r = 0, t = 0, since R fll/ R ,J,u = ^p7r 2 p 2 ~ T -4 and 
central proper time T ~ (— t) n . Notice that a does not approach unity as £ — > 00, but rather 
approaches a ~ 1.07. Likewise, mass grows linearly with distance, with m{r)/r — > 0.0596 
and Q — > 9.56 x 10~ 3 , so the solution is not asymptotically flat. An asymptotically-flat 



spacetime would require truncating the self- similarity at some large r [fLl j. Nonetheless, the 
self-similar solution is anticipated to represent the asymptotic behavior within 1Z (i.e., as 
r -»• 0, t -> 0). 

To confirm these expectations, we compute a series of models of spherically-symmetric 
collapse of radiation fluid, searching for a critical point and critical behavior. We adopt 
polar time slicing and radial gauge, retaining the form ([I]) of the line element, and assume 
a radiation fluid. We define W = aU l and U = aU r , with which the velocity normalization 



becomes W 2 = 1 + U 2 . Then, the equations of motion of the fluid, V V T^ V = 0, are 

0, (2) 
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and the pressure is p = |p. In these coordinates, the gravitational field equations C 
SnT r r and G t t = 87rT*i become 



(3) 



-tt + — (a 2 - l) = 4irra 2 p Uw 2 - 7 + 1 
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Anra 2 p "fU 2 + 7 — 1 



(4) 



(5) 



with 7 = 4/3 adopted in the models we discuss here. 

To construct a parameterized family of evolutions, we choose Cauchy data by specifying 
that the fluid be instantaneously at rest V = U /W = and have an energy density profile 



- 1 - -3/2 -2 / 2 / 2\ 

p = -7T 1 r r] exp(-r /r ), 



(6) 



parameterized by 77 and r . The total gravitational mass is M = ^Tq, so that r\ = 2M/r$ 
is a dimensionless measure of the strength of the initial gravitational field, while ro is the 
typical initial length scale. Equations (f|) and (H) are solved to complete the Cauchy data. 

The dynamical behavior of these models is consistent with expectations. For 7 = 4/3 
and this family of Cauchy data, the critical point is found in a bisecting search to be near 



Vc 



1.0188. For 7] > r/ c , the impending formation of a black hole is observed. As is 



well-known, polar time slices avoid penetrating regions of trapped surfaces [O] and asymp- 
totically approach the apparent horizon in collapsing spacetimes. When this occurs we have 
clear evidence of a black hole forming. Fig. 2 illustrates the form of the mass function 
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m(r) = \r{l — a~ 2 ) late in a supercritical evolution. The "kink" where r ~ 2m(r) locates 
the surface of the black hole and determines the hole's mass Mbh- Immediately beyond the 
kink, m(r) is nearly constant indicating an evacuated region that forms around the black 
hole. At larger radii, additional mass is evident, representing the fluid which did not collapse 
and is expanding outward. 

Black-hole mass Mbh is found to be well described by a power law (see Fig. 3) 

M BH ^ K( V - Vc f, (7) 

with a critical exponent (3 ~ 0.36. This critical exponent value is numerically indistin- 
guishable from those found in scalar wave collapse U and gravitational wave collapse [Q, 
doubtlessly reflecting a deep property of the gravitational field equations. 

The fluid is initially everywhere out of equilibrium. In near-critical models, the fluid in 
the inner region collapses while that in the outer region accelerates and expands outward. 
The collapsing central region is chased by a strong, ingoing rarefaction wave. As the radius 
R(t) of the inner region diminishes, the rarefaction wave causes the mass to diminish also, 
keeping m(R(t))/R(t) nearly constant. The radius R(t) of transition between collapse and 
expansion (edge of the rarefaction wave) decreases by several orders of magnitude in very- 
near-critical models. On fine scales, r < ro, the fluid and gravitational field assume a 
self-similar form (see Fig. 4). Furthermore, the self-similar form being approached is the 
locally self-similar solution outlined above. As Fig. 4 shows, ultimately the self-similarity 
of the collapse is broken if 77 ^ rj c . Only a precisely-critical model would approach the self- 
similar solution and retain this dependence as r — > 0, t — > 0. Near- but not precisely-critical 
spacetimes develop a self-similar region on scales r r$ but then lose self-similarity on a 
scale determined by proximity in parameter space to the critical point: r\ — K\rj — rj c f. 
Self-similarity will only be apparent if r% < r . 

Approach to self-similarity within a region spanning two disparate scales is a familiar 
concept in hydrodynamics termed "intermediate asymptotics" [|13|| . The cardinal example is 
the modified Sedov- Taylor blast wave problem. Besides the usual Sedov- Taylor parameters 
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of Eq, the blast wave energy, and po, the ambient gas density, the modified problem has a 
length scale r , within which the energy E is arbitrarily distributed initially, and has a small 
pressure po in the ambient gas. This small pressure introduces a second, new length scale: 
ri = (Eo/po) 1 / 3 3> r . These two additional dimensional parameters break the self-similarity 
of the blast wave on length scales r < r and r > r x but on intermediate scales the solution 
asymptotically approaches the Sedov- Taylor self-similar form. There is a direct analogue in 
the behavior of near-critical spacetimes on scales between r and r\. 

The self-similar solution determines the asymptotic properties of precisely-critical evolu- 
tions. Linear stability analysis should reveal the perturbative response of solutions nearby 
r/ c in parameter space and thereby perhaps provide an estimate of the critical exponent (5. 
The self-similar solution, and therefore any precisely-critical evolution, apparently contains 
a naked singularity. However, the precisely-critical evolutions are likely to be a set of mea- 
sure zero as there is extreme sensitivity to the initial conditions near r/ c . Furthermore, the 
dependence of the lapse along £ = oo (t = 0) is a ~ r l ~ l / n ^ r o.i29^ i nc ij ca ti n g the singularity 
to be shrouded in infinite redshift. 
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FIGURES 

FIG. 1. Self-similar radiation fluid collapse. Fluid and metric variables are plotted versus 
self-similarity coordinate £ = r/C{—t) n (with n ~ 1.1485): (1) radial metric function a(£) 
(solid), (2) velocity V(0 = aW/aU 1 (dotted), (3) = ^ P (dashed), and (4) lapse func- 

tion a = nrN(£)/t (dot-dashed). Note the collapsing interior and expanding exterior regions. 

FIG. 2. The mass function m{r) (solid) at a late time in a supercritical evolution. Polar time 
slices limit outside the horizon producing the kink near where the function approaches the curve 
(dotted) r = 2m(r), locating the surface of the black hole and determining the hole's mass. At late 
times, the fluid bifurcates and forms an evacuated region immediately beyond the hole, indicated 
by the decrease of Q = Anr 2 p (dashed) outside the hole. At larger radii, escaping fluid expands 
outward. 

FIG. 3. Critical behavior of black hole mass. For rj > 7] c , Mbh is determined and fit to a power 
law of the form (8). Here log of Mbh is plotted versus log of the critical separation r\ — r} c . The 
power-law fit is indicated by the solid line and is determined by i] c = 1.018828234, K = 3.27, and 
(5 = 0.36. We obtain black holes with masses down to Mbh — O.OOIMadm and the scaling law 
(8) is seen to hold over two orders of magnitude in mass. Simulations with discretization scales of 
3.3 x 10 -3 and 1.0 x 10~ 3 are included in the plot, using open and filled circles, respectively. 

FIG. 4. Intermediate asymptotic approach of a near-critical evolution to self-similarity. The 
metric function a(r) (light curves) at a series of times, equally spaced in logarithm of central 
proper time r = log(— T), obtained from a near-critical evolution. The rightmost curve is near 
the beginning of the simulation and resembles a{r) in the Cauchy data. As collapse proceeds, the 
peak of a(r) shifts to progressively smaller scales. The metric and fluid approach the self-similar 
solution depicted in Fig. 1. The self-similar form of a(£) (see Fig. 1) is overlaid at r = —0.5, clearly 
indicating approach to self-similarity. At later times, self-similarity is broken as a black hole begins 
to form and a — > oo on the horizon. 
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